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The well-known Bogoliubov expression for the spectrum of a weakly interacting dilute Bose gas 
becomes inadequate when the density or interactions strength are increased. The corrections to 
the spectrum due to stronger interactions were first considered by Beliaev (S. T. Beliaev, Sov. 
Phys.-JETP, 7:289, (1958)). We revisit Beliaev's theory and consider its application to a dilute gas 
with van der Waals interactions, where the scattering length may be tuned via a Fano-Feshbach 
resonance. We numerically evaluate Beliaev's expression for the excitation spectrum in the interme- 
diate momentum regime, and we also examine the consequences of the momentum dependence of 
the two-body scattering amplitude. These results are relevant to the interpretation of a recent Bragg 
spectroscopy experiment of a strongly interacting Bose gas. 



I. INTRODUCTION 

A dilute Bose Einstein condensate with sufficiently 
weak interactions is a many body system which can be 
described to a great degree of success by mean field the- 
ory. On the other hand, condensed-matter superfluids 
such as '^He are much more complex systems and diffi- 
cult to describe from first principles due to their high 
density and stronger interactions, leading to, in super- 
fluid Helium, a condensate fraction of only about 15% 
[l| . Superfluid Helium also exhibits the interesting prop- 
erty of a roton minimum in its excitation spectrum, which 
is absent in a weakly interacting dilute alkali gas. It is 
therefore of theoretical and experimental interest to try 
to bridge the two regimes by studying a Bose gas with 
stronger interactions. 

Dilute gases with stronger interactions can be stud- 
ied by considering perturbation theory corrections to the 
mean field theory. Lee, Huang and Yang (LHY), and 
Brueckner and Sawada 0, H, first determined the 
quantum depletion and the correction to the chemical 
potential of the ground state of a homogeneous dilute 
Bose gas. It was found that the leading corr ectio ns are 
proportional to the dimensionless parameter V na^ where 
n is the density and a the scattering length of the atoms. 
Beliaev Q first applied the methods of quantum field 
theory, specifically Green's functions, to the description 
of a system of bosons. We shall distinguish between 
the first order Beliaev's theory, which is equivalent to 
the Bogoliubov spectrum for the case of cont act i nterac- 
tions, and second order theory, which gives V na^ correc- 
tions to the excitation spectrum. These two levels corre- 
spond roughly to first and second order expansions of the 
Green's functions of the system in perturbation theory. 
A description of the quantum field theory of bosons and 
the Green's functions technique is available in textbooks 
0, Beliaev's work thus extends the LHY results for 
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the ground state to the excitation spectrum. 

A dilute Bose gas with tunable interactions can be 
realized experimentally with the utilization of a Fano- 
Feshbach resonance. A difficulty arises in that the gas be- 
comes more unstable with increasing interaction strength 
due to 3-body collisions. Nevertheless, in the shorter life- 
time of a more strongly interacting gas, it is still possible 
to observe relatively high frequency excitations (which re- 
quire short observation time, according to Heisenberg's 
uncertainty principle). Such an experiment has been re- 
cently performed in JILA [9]. 

Most studies of a dilute Bose gas assume an effective 
potential in the form of contact interactions. This model 
is usually very successful, but a correction is needed to 
describe excitations with higher momenta, due to the mo- 
mentum dependence of the scattering amplitude of two 
colliding atoms. In section II we shall examine how this 
correction can be incorporated in the first order theory 
of Beliaev Q . 

In section HI we revisit Beliave's second order theory, 
now assuming, for that purpose, a momentum indepen- 
dent scattering amplitude. Under this condition, Beliaev 
derived a closed form integral expression for the exci- 
tation spectrum. However, the integrals could only be 
evaluated analytically in the low and high momentum 
regimes. It is found that for low momenta (compared to 
1/^, where ^ is the healing length), the effect of stronger 
interactions is an upward shift in the excitation frequency 
(compared to the Bogoliubov spectrum), while at high 
momenta, the effect is to lower the excitation frequency. 
These results were also re-derived using the pseudopoten- 
tial method The interesting question that arises is: 
what happens in the intermediate regime? where is the 
transition point from upward shift to downward shift? 
This question is of particular relevance to the experi- 
ment [9|] which probed the spectrum at this intermediate 
regime. In this paper, we answer this question by eval- 
uating numerically the relevant expressions. To the best 
of our knowledge (perhaps somewhat surprisingly), this 
is the first time that Beliaev's theory's prediction for in- 
termediate momenta is examined. 
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II. EFFECT OF A MOMENTUM DEPENDENT 
SCATTERING AMPLITUDE 

Beliaev's original first order theory [H| already incorpo- 
rated, in general, the effect of a momentum dependent 
scattering amplitude. Here, we shall take a few more 
steps to specialize it to the description of a dilute gas 
with van der Waals interactions. We shall begin with 
some definitions. We use units h ^ 1 and also let the 
mass of the bosons be m = 1. 

Let n be the homogeneous total density of the gas, and 
no the density of the condensed part. In the original pa- 
per, Beliaev formulated his theory in terms of Uq. At the 
first order of theory, which we discuss in this section, the 
two are equal. Quantum depletion, giving rise to no < n, 
arise at the second order of the theory. Let ^p(r") be 
a wave-function solution of the Schrodinger equation of 
the two-body scattering problem with interaction poten- 
tial U (r), which behaves at infinity like a plane wave e*^ '' 
and an outgoing spherical wave. The scattering ampli- 
tude f{p',p) is related to the 5'(+)-function by 



fip',p)= J e-'P'-'^Uir)¥+\r)dr. 



(1) 



Notice that this definition differs by a numerical factor of 
—An from the usual one. The usual scattering amplitude 
is the value of f{p',p) at = |p| ('on the energy shell'). 
Eq. ([T]) thus generalizes the usual scattering amplitude to 
define an 'off the energy shell' amplitude. We shall also 
need a symmetrized amplitude defined as: 



(2) 



According to Beliaev's first order theory, the excitation 
spectrum of a dilute Bose gas is given by the dispersion 
relation: 



.(1) 
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where 6° = 



is the free particle kinetic energy. The 

superscript (1) in ep'' indicates first order. In the case 
of a gas with van der Waals interactions, it is possible 
to derive analytic expressions for the scattering ampli- 
tude. At low energies, s-wave scattering is the domi- 
nant process so that for on shell scattering p = p' we 
may define an isotropic amplitude /(p) = f{p',p)- From 
the usual definition of the scattering length a we have, 
/(O) = /o = 4™. 

The scattering amplitude is related to the s-wave phase 
shift 5q[p) as 



p cot So [p) — i 



(4) 



From the standard 
pcotSoip) = + ^"^eP 

fective range. Thus 

Hp) = ^ 



effective-range 
^ where 



expansion, 
is the ef- 



^arep"^ 



lap 



(5) 



For neutral atoms with van der Waals interactions, the 
following useful relationship exists between the scattering 
length g, the effective range re, and the Ce coefficient 
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1 + [1 - Xe(a//?6)]' 



(6) 



where /Jg = {mCe/h^)^^''' is a length scale associated with 
the van der Waals interaction, and Xe = [r(l/4)]^/(27r), 
with F being the usual Gamma function. For clarity, we 
have written here explicitly factors of h and m. Typical 
values of Pe for alkali atoms are about 100 bohr. 

By extending the ideas of [ll| , we find that the off-shell 
amplitude f{p,0) which appears in Eq. ([3]), is given, for 
a van der Waals interaction, by: 



1. 



fip,0)^Ana{l + -bp' + ...) 



with 



^--30^' 



5\/2/36/c 



, r(9/4) \ 
r(7/4); 



(7) 



(8) 



We have checked this analytic expression for the off shell 
amplitude with numerical calculations using a model po- 
tential of *^Rb, with good agreement. 

A few notes are in place. First, the appearance of 
the factor i in the expression §) is due to mov- 

ing to the center of mass system of an excited particle 
with momentum p and a condensate particle with zero 
momentum. Second, in using Eq. ([5]) with Eq. ([3]) one 
should note that the scattering amplitude has an imag- 
inary part. This gives rise to an imaginary part of the 
energy describing correctly, for high momentum, the de- 
cay of a particle excitation due to the scattering of the 
particle off the rest of the condensate's atoms. However 
at low momenta (the quasi-particle phonon regime), the 
decay is actually suppressed, and is only seen in the next 
level of theory (Beliaev's second order theory). At first 
order of the theory it is sufficient to take the real part of 
the scattering amplitude when using it to calculate the 
spectrum. The off shell scattering amplitude f{p, 0) is 
always real. 

for pa <^ 1, eq. ([3]) reduces to the familiar Bogoliubov 
spectrum: 



2no/oe^. 



(9) 



By increasing pa one enters the regime where finite- 
momentum corrections are necessary. For a fixed p, this 
could happen due to increasing a. If the scattering length 
is larger than the van der Waals length scale, i.e, a^ f3e, 
then according to Eq. ([6]), re « 1.4/36, while according to 
Eq. JHI), 6 « 0.44/3|. If the condition p/3e < 1 is stiff sat- 
isfied (which is normally the case), then the correction 
to the off shell amplitude /(p, 0) is very small and one 
can take it to be equal to a. Furthermore, in Eq. the 
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FIG. 1: Change in line shift due to the momentum depen- 
dence of the scattering amplitude. The three curves corre- 
spond to different values of the parameter a/^ = vovmc?. 
The values are (solid black), 0.2 (dotted blue), and 0.4 
(dashed green) 



with the excitation spectrum rather than the ground 
state, and includes other second order processes such as 
Beliaev's damping (the decay of a quasi-particle into two 
lower energy quasi-particles). Beliaev summed the rele- 
vant Feynman diagrams which contribute to the Green's 
function of the system up to second order. The excita- 
tion spectrum is then found by the poles of the Green's 
function in the {p,p'^) space, with p momentum and p'^ 
energy. Mohling and Sirlin [lO| arrived at the same exci- 
tation spectrum using a more straightforward application 
of the LHY pseudopotential method. We shall briefly 
state here the results in the language of Green's func- 
tions 0, 3 ■ I'^ this section we use the 4-momentum no- 
tation p = (p;p°). In a Bose condensate there are two 
Green's functions, a regular Green's function Gii(p), and 
an anomalous Green's function Gi2{p). Correspondingly, 
there are two self-energies, a regular self-energy I]ii(p) 
and an anomalous self-energy Si2(p). Also, we have the 
free particle Green's function G^{p). The Dyson's equa- 
tions of this system can written in the matrix form: 



effective range contribution due to non-zero is smaller 
by a factor p/3g than the universal universal factor ipa. 
Therefore, we arrive in a final universal momentum de- 
pendent dispersion relation, valid for p/Jg ^ 1: 



G{p) = G\p) + G\p)i:{p)G{p), 




1 + {pay/ 4 



-n2(47ra)2.(10) 



At this level of the theory, we may replace uq in Eq. ([TO] 
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with n. We define an energy shift E — e 
the difference between the excitation energy of a quasi- 
particle with momentum p in the condensate, and the 
kinetic energy of a free particle with the same momen- 
tum. In Fig. ([1]) we plot the energy shift as a function of 
the dimensionless parameter where ^ — is the 

healing length. The shi ft is sh own for different values of 
the parameter a/^ = VSTma^. The curve corresponding 
to a/^ — should be interpreted as taking the limit of 
V na^ — > (notice also that the vertical axis is propor- 
tional to na), and gives the usual Bogoliubov spectrum, 
with the asymptotic shift value of Aima for large momen- 
tum. Increasing a/^ gives rise to observable deviations 
at high momentum. At small momenta the shift is linear 
in momentum and the slope gives the speed of sound in 
the condensate. Since the change in the real part of the 
scattering amplitude is proportional to p^, the slope of 
all curves is the same at small momenta, and the speed 
of sound is unchanged by this effect. 



III. EFFECT OF QUANTUM DEPLETION 

The second order of Beliaev's theory is associated with 
the LHY quantum depletion, but is concerned mainly 



with 



G{p) 



G\p) 



Gii(p) Gi2(p) 
Gi2(-p) Gii(-p) 



GO(p) 
GO(-p) 



T.ii{p) Si2(p) 
T.i2(-p) Sii(-p) 



(11) 



(12) 



(13) 
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The self-energies can be calculated in perturbation the- 
ory using Feynman's diagrams. Then, the Green's func- 
tions can be obtained by solving the algebraic matrix 
equation (fTTjl. and their poles give the excitation spec- 
trum (both kinds of Green's functions have the same 
poles). Following the original papers, we assume in this 
section (in contrast to the previous one) a momentum 
independent scattering amplitude, valid for pa 1, i.e, 
/o ~ Attu. The first order of theory gives: 



Si2(p)(^) - nolo. 



(15) 
(16) 



When these expressions are used to calculate the excita- 
tion spectrum via the poles of the corresponding Green's 
functions, one retrieves the ordinary Bogoliubov spec- 
trum, Eq. (HI). 

The second order contributions (additive to the first 

(2) (2) 

order), written T,\^ and '£,[2 , are: 
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dq 



g(i)g(i) 



Q-{q,k) 



,(1) ,(1) 



where k = p q, and the functions i?, are: 

i?(q,fe) = 2e°4-26WeW+n^/„^ 



Q^{q,k) 



.(1) 



i<5 



Q+{q.k) 



i(5 



I — 



f(i).o _ (1) 



(17) 



(18) 



The id factors in Eq. (fT7|) are needed for convergence and 
should be understood in the usual sense of the limit 5 
0"*". By solving the Dyson's matrix equation with these 
self energies, the Green function Gn can be expressed in 
the convenient form (showing its poles): 



Gii(p) = 



A,, 



Br, 



A-{p) pO + ei''^ +A+{p) 



, (19) 



where Ap and Bp are certain functions of p (with no 
poles), and A^(p) are second-order corrections 



A^(p) = -|^(Sii(p) + Sii(-p)-2^) 



(2) 



2e 



nofo 



(20) 



^(Sn(p) + En(-P) -2^l- 2^^)^''^ ± 



2e« 

±i(En(p)-Sn(-p))(^), 

with /i^^) the second order correction to the chemical po- 
tential, 



^(2) ^ (5/37r^)A/no/3no/o 



(21) 



Inspecting the poles of the Green's function, the exci- 
tation spectrum in this approximation is then given by: 



6p = eW+A-(p;e«). 



In general, there is no known analytic solution to the 
integrals in Eq. (fT7|) . Beliaev derived an analytic solution 



in the the limit of small momenta <C 1: 

1 



Py nfo[l + 

(pf < 1). 
6407rn ' 



(23) 



This result is expressed here in terms of the total density 
n rather than the condensate density uq. The relation 
between n and no is: 



- no = no J no fg/ {Sir ). 



(24) 



According to Eq. ((23)) . the second order effects give rise 
to a upward shift in the excitation energy with linear 
dependence in momentum, and thus a shift in the velocity 
of sound: 



c = y fon 



1 



1 + 



(25) 



This result for the speed of sound, derived above from 
the slope of the one-particle excitation spectrum at low 
momentum, is identical to that obtained by considering 
only the LHY expression for the ground state energy: 



V 2m 



16 
15^ 



(26) 



(22) where we explicitly included factors of m and h, and using 
the thermodynamic relations P = — (f^)^ for the pres- 
sure P, and mc^ = §^ f^- Thus, from the macroscopic 
thermodynamic point of view, the increase in speed of 
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sound is traced to the equation of state becoming stifFer 
due to quantum depletion. Eq. also has an imagi- 
nary part which is due to the decay of one quasi-particle 
into two other quasi-particles. 

On the other hand, Mohling and Sirlin give an 
expression for the shift at high momenta, with ^ 1, 
(but, at the same time, pa <^ 1). The result may be 
written in the form: 

Cp = el + 2nfo- fiLHY - (27) 

- tna^ (K»l), 

with the LHY chemical potential : 

t^LHY = nfo + ^ , (28) 

and the cross section for two identical bosons 

cr = Stto^. (29) 

We may understand the result of Eq. ([27]) as follows: 
the excitation energy is the difference between final state 
and initial state energies. The initial state energy of the 
(quasi-)particle is /x, and the final state energy at high 
momentum is the free kinetic energy plus the sum of 
direct and exchange interactions of the excited particle 
with all other, low momentum particles, 2no/o- In ad- 
dition, the excitation decays due to collision of the high 
momentum particle with all other particles, with cross 
section a and center of mass momentum p/2. Note, that 
Eq. ((27l) only considers the high momentum {p^ ^ 1) 
effect due the many body interactions. The correction 
due to the momentum-dependent scattering amplitude, 
as described by Eq. (fTU|) can be neglected here as long as 
pa <C 1. The required conditions, 1/^ <Cp<C 1/a, can 
in principle be satisfied simultaneously for a sufficiently 
dilute gas. 

What happens in the intermediate range, « 1 ? 
We have calculated numerically the excitation energy 
Eq. (P^ . This was done by deforming the path of inte- 
gration from the real line into the complex plane around 
the singularities indicated by the i6 factors. The re- 
sult is shown in Fig. The plot shows the quantity 

3?(ep — ep ^), i.e the difference between the real part of Be- 
liaev's excitation energy and the Bogoliubov energy. At 
low momenta there is an upward linear shift consistent 
with Eq. . This linear shift is consistent with the cor- 
rection to the speed of sound derived by LHY theory. The 
linear behavior holds up to ~ 0.6. At high momen- 
tum there is a down shift with very slow 1/p convergence 

towards the asymptotic value ^fon\J /qh/'k'^ consistent 
with Eq. (j27p . In between, we find the transition point 
at p « 6.18^"^, where the second order contribution is 
zero, so that at this point the excitation energy to second 
order is identical with the first order Bogoliubov's theory. 




FIG. 2: Solid black line: change in line shift due to quantum 
depletion and other second-order processes. The correction 
is additive to the solid black line in Fig. ([l]). The amount of 
correction is proportional to 8^/naFjTv, as indicated on the 
vertical axis. Blue dashed line: the linear asymptote at small 
momenta which is consistent with LHY theory. Red dotted 
line: the asymptotic value at high momenta. 
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FIG. 3: Cross section for quasi-particle decay. 

The imaginary part of the excitation energy gives rise 
to quasi-particle damping with an effective cross-section 
CTp = — 29(^) ranging from at p = to Stto^ a.tp£_ ^ 1. 
It is shown in Fig. This damping cross-section has 
been calculated before using another method and mea- 
sured in experiment [l^ . The two methods of calculation 
give identical results. 



A. Discussion and Conclusions 

In this paper we have revisited Beliaev's theory with 
emphasis on two mechanisms for a shift in the excitation 
spectrum of a Bose gas from the usual Bogoliubov the- 
ory. One is due to the momentum dependence of scatter- 
ing amplitude, when the condition fca <C 1 is not valid, 
which generally causes a downward shift in the excita- 
tion frequency at higher momentum, as demonstrated in 
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Fig. ([T|). The other is a purely many body effect, occur- 
ring when the condition voTmo^ ^ 1 is not obeyed. This 
effect is associated with quantum depletion due to inter- 
actions as well as with quasi-particle decay. As shown in 
Fig. ([2). It gives rise to upward shift at low momentum 
and downward shift at high momentum. 

We note that the two separate effects were calculated 
under exclusive conditions. The effect due to momentum 
dependent scattering amplitude {pa « 1) was calculated 
in first order of the interactions, valid for \/ Snna^ ^ 1 . 



On the other, the second order effects describing VSnna-^ 
corrections assumed pa <ti 1. Currently, there is no 
known result that holds in the regime where both the 
momentum dependence of the scattering amplitude and 
quantum depletion due to strong interactions are impor- 
tant. The physical interpretation that follows Eq. ([27|) 
suggests that at high enough momentum the shift due to 
these two effects should be additive and combine to give 



asymptotic value. The Bragg spectroscopy experiment in 
Ref. Q explored the regime where both effects discussed 
above are important, and showed evidence for a down- 
ward shift in the excitation energy of a strongly inter- 
acting gas. But for the measurements with the strongest 
interactions, which showed the largest shifts, the exper- 
iment had p^ ~ 2, where Fig. ([2|) shows upward rather 
than downward shift. On the other hand, the prediction 
of down shift due to the 2-body momentum dependence 
of the scattering length seems to go some way towards 
explaining the observations. 

Finally, we note that the dynamic structure factor is 
a quantity more directly related to the actual measure- 
ment process of a Bragg experiment QjllBl ■ ^ partial at- 
tempt at incorporating second order effects due to strong 
interactions in the dynamic structure factor is found in 
Ref. d^. It would be of interest to complete this pro- 
gram. 



el + 2n^ 



\P\ 



fs{p/2) 

(K » 1) 



^J-LHY 



(30) 
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